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We study a generalized multichannel single-impurity Kondo model, in which the impurity spin is 
described by a representation of the SU(iV) group which combines bosonic and fermionic degrees 
of freedom. The impurity spin states are described by Abrikosov pseudofermions, and we make 
use of a method initiated by Popov and Fedotov which allows a proper handling of the fermionic 
constraint. The partition function is derived within a path integral approach. We use renormal- 
ization group techniques to calculate the /3 scaling function perturbatively in powers of the Kondo 
coupling constant, which is justified in the weak coupling limit. The truncated expansion is valid 
in the overscreened (Nozieres-Blandin) regime, for an arbitrary SU(iV) group and any value of the 
parameters characterizing the impurity spin representation. The intermediate coupling fixed point 
is identified. We derive the temperature dependence of various physical quantities at low T, con- 
trolled by a unique critical exponent, and show that the physics of the system in the overscreened 
regime governed by the intermediate coupling fixed point is characterized by a non-Fermi liquid 
behavior. Our results are in accordance with those obtained by other methods, as Bethe ansatz 
and boundary conformal field theory, in the case of various impurity spin symmetries. We establish 
in a unified way that the Kondo models in which the impurity spin is described successively by 
a fundamental, symmetric, antisymmetric and mixed symmetry representation yield all the same 
low-energy physics in the overscreened regime. Possible generalizations of the analysis we present 
to the case of arbitrary impurity spin representations of SU(iV) are also discussed. 



I. INTRODUCTION 

In spite of its twenty-five years of history, the mul- 
tichannel Kondo modeli*^ is still a subject of current 
interest in several fields of study of condensed matter 
physics. Indeed, this simple impurity model exhibits 
a non-Fermi liquid (NFL) behavior. Originally intro- 
duced to describe metallic systems containing dilute mag- 
netic impurities, its applications have been consider- 
ably extended. For instance, it is related to two-level 
systemsi 3 i 4 i 5 It has been proposed^ that materials such 
as Pbi-^Ge^Te and Ki^Li^d could be modeled using 
a multichannel version of the Coqblin-Schrieffer models 
i.e., a SU(iV)xSU(A') Kondo model. Here N indicates 
the number of spin degrees of freedom, and K is the 
number of channels for the conduction electrons. The 
multichannel Kondo model is also believed to be relevant 
to the problem of mesoscopic quantum dots^S Several 
additional possible realizations concern the context of 
heavy- fermion compounds. Recent experiments in such 
systems have shown the existence of a quantum phase 
transitionifi from a magnetically disordered to a long- 
range magnetic ordered phase, driven by a change in 
chemical composition, pressure or magnetic field. One 
can cite the alloy CeCu5.9Auo.1ri the behavior of which 
in the disordered phase close to the quantum critical 
point is clearly of NFL-type: the specific heat C depends 
on T as C/T ~ — ln(T/T ), the magnetic susceptibility 
as x ~ l~aVr, and the T-depcndcnt part of the resistiv- 



ity as Ap ~ T instead of C/T ~ Const., \ ~ Const, and 
A/9 ~ T 2 as predicted by the Fermi liquid theory^ Simi- 
lar types of behavior have been observed in other systems 
such as CePd2Si2r 3 - and YbRh 2 Si2r^ among others (for 
an extensive survey of the experimental situation con- 
cerning heavy-fermion compounds, the reader is refcred 
to the review article of Steward). The experimentally 
observed breakdown of the Fermi liquid theory addresses 
fundamental questions about the possible formation of 
novel electronic states of matter with new types of ele- 
mentary excitations resulting from the presence of strong 
correlations among electrons. 

Since it has been proposed the multichannel Kondo 
model has been studied intensively by a variety of the- 
oretical methods, such as the numerical renormaliza- 
tion group (RG)r^ Bethe ansatz^S*^ conformal field 
theory; 20 i 21 i 22 i 23 and perturbative RGi^ Detailed dis- 
cissions concerning the various experimental applica- 
tions and theoretical aspects of the multichannel Kondo 
model can be found in the book of Hcwsonr& and the 
review article of Cox and Zawadowski^i 

Recent experiments performed in the heavy-fermion 
compound CeCus.gAuo.i by inelastic neutron scattering 
have shown the existence of an anomalous uj/T scal- 
ing law for the dynamical spin suceptibility at the an- 
tiferromagnetic wavevector which persists over the en- 
tire Brillouin zoneiSLSi Such a result indicates that 
the spin dynamics are critical not only at large length 
scales, as studied in the traditional itinerant magnetism 
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picture f 3 ^ 3 ^* 3 ^ 3 ^ 3 * but also at atomic length scales. It 
strongly suggests the presence of locally critical modes 
beyond the standard spin-fluctuation theory. An ap- 
proach which has been put forward to describe such a 
local quantum critical point is based on extended dy- 
namical mean field theory, with encouraging results like 
the prediction of a scaling law for the dynamical spin 
suceptibility^ Another theory is based on the general- 
ization of the spin operators symmetry, from SU(2) to 
SU(N), to describe the Kondo impurity spin in terms 
of a mixed fermionic-bosonic representation^ The idea 
behind this mixed formulation is to account for both the 
extended quasiparticle excitations and the locally criti- 
cal modes using the fermionic and the bosonic part of the 
spin, respectively. 

In order to fix ideas, let us start with the antifcrro- 
magnctic single-channel SU(2) Kondo impurity model. 
It is well known that within a RG analysis j2LSi2M2i the 
flow takes the Kondo coupling Jk all the way to strong 
coupling (i.e., to infinity). An important aspect in the 
discussion of the breakdown of the Fermi liquid theory 
is related to the question of the stability of the strong 
coupling (SC) fixed point. In this case the system flows 
to a SC fixed point which is stable and the associated be- 
havior of the system is that of a local Fermi liquid^ 1 ^ 2 *^ 3 - 
The situation is rather different when one considers sev- 
eral channels of conduction electrons. In the case of a 
spin S interacting antiferromagnetically with the spin of 
conduction electrons belonging to K different channels, 
Nozicrcs and Blandini have shown that the multichannel 
Kondo model can lead to two very different situations de- 
pending on how K compares to 2S. Their calculation is 
based on a second order perturbation theory in the hop- 
ping amplitude of the conduction electrons, around the 
SC fixed point. They analyze their results by deriving 
an effective coupling between the spin of the composite 
formed by the impurity dressed by the conduction elec- 
trons in the SC limit, and the spin of the conduction elec- 
trons on the neighboring sites. On the one hand in the 
undcrscreened regime, when K < 2S, the effective cou- 
pling is found to be ferromagnetic and the SC fixed point 
is stable. In this regime an effective spin, (S — K/2), is 
formed, resulting from the screening of the impurity spin 
by the conduction electrons located on the same site. The 
system, described by the SC fixed point, behaves as a lo- 
cal Fermi liquid. On the other hand in the overscreened 
regime, when K > 2S, the effective coupling is found 
to be antiferromagnetic and hence the SC fixed point is 
unstable. The instability of the SC fixed point in this 
regime constitutes a hint for the existence of an interme- 
diate coupling (IC) fixed point, characterized by a finite 
critical value of the Kondo coupling and leading to a NFL 
excitation spectrum. 

In a previous work^i^ using a second order pertur- 
bation theory similar to the analysis of Nozieres and 
Blandin, we have considered an antiferromagnetic single- 
channel SU(iV) Kondo model involving a generalized im- 
purity spin, realized by a combination of 2S bosonic and 



q fermionic degrees of freedom. Preceding studies* 6 - have 
shown that the SC fixed point associated with this model 
includes two regimes, one of them characterized by an 
instability for a certain class of mixed symmetry impu- 
rity spins. More precisely in the large- TV limit the SC 
fixed point is stable and describes the low-energy physics 
of the one conduction channel model when q < N/2, but 
becomes unstable for q > N/2. This result suggests that 
the instability of the SC fixed point can arise from the 
symmetry properties of the impurity spin, without any 
introduction of additional degrees of freedom like con- 
duction electron channels, coupling anisotropy or else. 
We have investigated how the system behaves, in partic- 
ular when the number of conduction electrons (filling) on 
the neighboring site of the impurity located at the ori- 
gin varies. We have expressed our results in terms of an 
effective hamiltonian describing the interactions between 
the impurity dressed by the conduction electrons at the 
origin (as obtained in the SC limit), and the conduction 
electrons located on the neighbouring site. By deter- 
mining the effective couplings, we have shown that the 
SC instability, resulting from the change in sign of the 
effective spin interaction, is related to the change from 
repulsive to attractive of the effective charge interaction. 

In the present work we consider the multichannel ver- 
sion of the model described above, i.e., a SU(AT)xSU(.K') 
Kondo model, in order to achieve an even richer situation. 
In the overscreened regime, the magnetic degrees of free- 
dom of the impurity arc completely quenched, and the SC 
fixed point is unstablei We present a method allowing 
the identification and the characterization of the IC fixed 
point in the overscreened regime. Our study is based on 
a perturbative RG method pioneered by Anderson, the 
"poor man's scaling" approach ; 26 ' 37 Its basic principle 
is to reduce the conduction electron band width, D, by 
integrating out the high-energy intermediate states, and 
letting the low-energy physics properties unchanged. The 
main difficulty for the description of spin systems is that 
spin operators obey neither Fermi nor Bose statistics, 
leading to the absence of a Wick's decomposition theo- 
rem in this context. Consequently, conventional pertur- 
bative methods such as a Feynman's diagram approach 
are prohibited. The Gaudin theorem 4 ^ existing for SU(2) 
generators does not allow the construction of simple di- 
agrammatic techniques directly for spin operators. We 
tackle this problem by making use of a method initiated 
by Popov and Fedotov^S In the framework of a fermionic 
description of the impurity spin states, this powerful ap- 
proach allows a proper handling of the fermionic con- 
straint. It introduces a set of auxiliary chemical poten- 
tials, which eliminates the unphysical Hilbert sub-space 
without removing the true spin states. The overscreened 
regime of the Kondo model has been studied intensively 
by several methods, for various impurity spin representa- 
tions in both the SU(2) 24 i 25 i 5Q i 51 i 52 and SU(AQ 9 i 53 i 54 i 55 i 56 
spin groups. Our work considers several of these impu- 
rity spin representations as particular limits of the mixed 
symmetry formulation, enabling direct comparisons with 
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the results obtained by other methods. 

This paper is organized as follows. In Section HT1 the 
model is introduced, and we present the formalism lead- 
ing to the derivation of the partition function within a 
path integral approach. We adopt the Abrikosov pseud- 
ofermion formulation to describe the impurity spin states. 
The handling of the fermionic constraint, performed by 
making use of the Popov-Fcdotov method, is discussed. 
In Section IIIII we use a diagrammatic method to de- 
rive the conduction electron self-energy perturbatively in 
powers of the Kondo coupling constant, i.e., in the weak 
coupling limit. The imaginary part of the conduction 
electron self-energy is proportional to the electronic scat- 
tering rate, which is invariant under RG transformation 
as any physical quantity. This invariance is exploited in 
Section lTVI to compute the /3 scaling function by making 
use of perturbativc RG arguments, including the Callan- 
Symanzik equation. This leads to the identification of the 
IC fixed point in the overscreened regime. The temper- 
ature dependence of various physical quantities at low 
T, controlled by a unique critical exponent, is derived. 
The NFL nature of the physics governed by the IC fixed 
point is also discussed, and our results are compared with 
those of previous works. Section gives a summary of 
the method and main results, and is devoted to conclu- 
sions. 

The appendices contain the technical details of the 
calculations. In Appendix ^ we give some properties 
followed by the SU(iV) group generators. Appendix [B] 
briefly explains the principle of the method pioneered by 
Popov and Fedotov. This method requires a particular 
derivation of the conduction electron and impurity spin 
pscudofcrmion free propagators, which is performed in 
Appendix [Cl Appendix IdI contains the calculation of the 
first order term spin factors showing up in the conduc- 
tion electron self-energy expansion. In Appendix [E] we 
show exhaustively how to compute the contribution to 
the electronic self-energy for a given example of diagram. 
The other diagrams can then be calculated in a similar 
way. 



II. GENERALIZED KONDO MODEL AND 
FORMALISM 



In this Section we present the model, and the methods 
used to derive the partition function which is obtained in 
the functional integration formalism. In the language of 
Young tableaus, the impurity spin is described by a L- 
shaped representation, as shown in Fig. ^ The bosonic 
degrees of freedom correspond to the horizontal line (first 
row) containing 2S boxes, while the fermionic ones are 
included in the q boxes of the vertical line (first column). 
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FIG. 1: Young Tableau description of the L-shaped irre- 
ducible representation {T} associated with the mixed sym- 
metry impurity spin, realized by a combination of bosonic 
and fermionic degrees of freedom. 

A. Generalized multichannel SU(iV) 
single-impurity Kondo model 

The Hamiltonian describing the Kondo model we con- 
sider is 

H K =H + V Kl (1) 

where Ti.o is the free electronic part and Vk gives the in- 
teraction between the generalized impurity spin and the 
conduction electrons. A SU(Y) (for N > 2) spin degener- 
acy is required for both the impurity and the conduction 
electron spins, and we suppose that there exists K (> 1) 
electronic conduction channel(s) 

N K 

H o = J2J2J2 e k cl aX c k ,^ x , (2) 

k cr=lA=l 

where c\, x is the creation operator of a conduction elec- 
tron with momentum fe and energy e k , <J and A are the 
spin and channel indices, respectively. The second term 

V K = J K Si ■ S e (0) , (3) 

corresponds to the antifcrromagnctic Kondo interaction 
Jk (> 0) between the impurity spin Si located at the 
origin (site 0) and the spin of the conduction electrons 
S e (Q) located at the same siteS 

The impurity spin states corresponding to the irre- 
ducible representation {T} can be expressed in terms of 
Abrikosov pseudofermions^ We denote the dimension of 
{r} by d{r}- Its general expression is given in Appendix 
lAl Each impurity spin state is described by a fermion. 
Therefore we have species of pseudofermions located 
at the origin which are respectively created and destroyed 
by the following operators 

4 , / 7 , 1 < 7 < d {r} . (4) 

The number of pseudofermions at the impurity site is 
subject to the standard constraint 

d lr y 

"/ = E 4/7 = 1 ■ ( 5 ) 

7=1 
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FIG. 2: Young Tableau description of the fundamental (irre- 
ducible) representation {a} associated with each conduction 
electron. 



According to the structure of the SU(iV) Lie group^ the 
impurity spin operator owns (iV 2 — 1) components 



ka<n 2 -i 



with 



7,7'=! 



(0) 



where T A corresponds to the (A^ 2 — 1) generators of the 
SU(iV) spin group. A crucial technical aspect has to be 
pointed out here. Due to the shape of the generalized 
impurity spin that we consider which is not in the fun- 
damental SU(iV) representation (i.e., a single box repre- 
sentation in the language of Young tableaus), and to the 
Abrikosov pseudofermion formalism that we adopt, the 
{T A } 1<A<N 2_ 1 S\J(N) generators are not expanded as 

(N x N) matrices as usual but as (rf{r} x <^{r}) ones. Such 
a situation already occurs in the SU(2) group if instead 
of considering a spin 1/2 corresponding to the SU(2) fun- 
damental representation, we work with a spin 1, 3/2, ..^H 
The properties of the \T A } matrices which we will need 
are given in Appendix IAl It is to be noted that this way 
of writing the impurity spin states, although equivalent, 
is completely different from the method that we used in 
our previous work4^*^ In the latter case, the impurity 
spin states were obtained as a combination of (25 + g — 1) 
particles, each box of the L-shapcd representation being 
described by either one boson or one fermion, and the 
SU(iV) generators we had to consider corresponded to 
the fundamental representation. However the approach 
we now consider in terms of pure pseudofermions is much 
more appropriate to the functional integration formalism 
we develop. 

As shown in Fig. |3 each conduction electron is de- 
scribed by the fundamental representation of SU(iV). 
The spin of the conduction electrons at the origin jf? e (0) 
is given by 



S e (0)= S A (0) 



1<A<N 2 -1 



with the (A^ 2 — 1) spin components given by 

N K 



'=1 A=l 



where cr A corresponds to the (N 2 — 1) generators of the 
SU(AT) spin group. In the previous expression c' , A (0) 



is the creation operator of a conduction electron of spin 
color a in the channel A at the impurity site, which can 
be rewritten in the momentum space as 



(7) 



where M s is the number of sites of the lattice. Using Eq. 
10, the different components of the conduction electron 
spin can be expressed as 



1 



N K 



<r A , a c k ^ x . (8) 



S k,k' er.er'=l A=l 



The SU(7V) generators for the conduction electron spin 
correspond to the fundamental representation, there- 
fore the corresponding tracclcss hcrmitian matrices 
{cr A } 1<j4<Ar2 l have (N x N) components and satisfy 
the following normalization equation 

T*[ct a <t b } = 1 -6 ab . 

In the well known SU(2) and SU(3) groups, a A = t a /2 
where {t" 4 } are the Pauli and Gell-Mann matrices, re- 
spectively. 

By making use of Eqs. J2J , © , © and © we ob- 
tain the Hamiltonian related to the multichannel SU(7V) 
x SU(if) single-impurity Kondo model for a L-shapcd 
Young tableau representation of the impurity spin 



N K 



E E E £ fe4, CT ,A c fe^, 



k er=lA=l 
N 2 -l 



Jk 
AC 



N K 



E (E E E c Iv,a<. c *^ 

A=l k.k' <7,<t' = 1 a=1 



( t A \ 



(9) 



7,7 =1 



B. Functional integration method - Handling of the 
fermionic constraint 



The partition function corresponding to TIk @ is 



Z = Tr spin exp [ - (3H K ] 



(10) 



with (3 = 1/fceT. The difficulty in computing Z comes 
from the constraint expressed in Eq. J5J which requires to 
perform the trace Tr sp i n over the physical states involv- 
ing one pseudofermion at the impurity site only. However 
the Pauli principle allows the number of pseudofermions 
at the origin to vary between and <i{r}- More explic- 
itly, the representation of the impurity spin with pseud- 
ofermion operators enlarges the dimensionality of the 
Hilbcrt space in which these fields arc acting, compared 
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with the physical space appropriate for the true spin op- 
erators. The unphysical states displayed in the pseud- 
ofcrmion formulation arc then eliminated with the aid 
of the constraint. Therefore the evaluation of the trace 
in Eq. 1|10|) can not be simply performed using standard 
perturbation approaches for instance but requires the use 
of further techniques. 

In order to settle this problem, we adopt a method 
pioneered by Popov and Fedotov^ As explained in Ap- 
pendix^ one can introduce a set of auxiliary (imaginary) 
chemical potentials {i^ n }i<ri<d {r} governed by a discrete 
distribution function P{n} 



1 "{r} 

P{n} = - — Y] exp [i/3/x n ] 6(fi + i/x«) , 



as 



l3d {T} 



(2n-l) 



(11) 



A crucial property arising from Eq. Hlf) is that the imag- 
inary chemical potentials depend on temperature. The 
partition function can be rewritten as 

Z = Jvn P{/i} Tr Fock exp[ - p(H K - fihf)] , (12) 

where the trace Trp ck is performed over the entire Fock 
space corresponding to the fermionic impurity spin oper- 
ators without any constraint, namely including any num- 
ber of pseudofermion(s) allowed by the Pauli principle. 
The contributions of the unphysical states (rif ^ 1) to 
the partition function automatically cancel out one with 
other by making use of the imaginary chemical poten- 
tials. With Eqs. JTTJ and we get 



1 a {r} 

Z = - ^2 ex P [i/^™ 



l {r} — ! 



(13) 



xTr Fock exp 



/3(H K + ifJ-nfif) 



In the latter equation, in contrast with Eq. (|l(Jfl of Z, the 
trace showing up in the calculation has to be performed 
in the grand canonical ensemble, which allows to use 
standard techniques, namely the functional integration, 
Green's function method and diagrammatic expansion. 
In a path integral formalism at finite temperature jSLSS 
each Fock trace is given by 



Tr Foc kexp 

/ Vc*VcVf*Vfexp 



(3(Hk + ifJ-nfif) 




(14) 



dr{£„(r)+W K (T)} 



where c*, c, /*, / are Grassmann variables associated 
with the fermionic operators c\ c, f\ /, respectively, 
and C n (r) is the Lagrangian of free particles expressed 



N K 



£«(t) = Y Y 4,<r,A( T )^ c fc^A(r) 

k tr=l A=l 
rf{r } 

+ ^/ 7 *(r)(a T + i Mn )/ 7 (r) . (15) 

7=1 

One can notice from Eq. H15|) that a direct consequence of 
the introduction of the auxiliary chemical potentials is a 
shift of the Matsubara frequencies related to the pscudo- 
fermions by a value of (— i/i„). 

Using Eqs. (|14f> and (|15fl . we finally obtain the expres- 
sion of the partition function 



Z = V exp [i(3fi n ] [ Vc*VcDf*Vf (16) 

d m ±-J J 

/ dT{C n (T)+H K (T)} 

Jo 



III. 



x exp 



CALCULATION OF THE CONDUCTION 
ELECTRON SELF-ENERGY 



In this Section, we expand the conduction electron self- 
energy perturbatively in powers of the Kondo coupling 
constant Jk- This means that we consider the weak cou- 
pling limit. We closely follow the method developed by 
Gan et alM& in the special case of the SU(2) x SU(X) 
Kondo model, involving an impurity of spin 1/2. We use 
diagrammatic techniques SSa 



A. Construction of the electronic self-energy 
diagrams 

Here we summarize the rules allowing to derive the 
self-energy of the conduction electron, using a diagram- 
matic pcrturbativc expansion of the partition function in 
Jk- A detailed method to build the diagrams associated 
with the electronic self-energy in the framework of the 
Kondo model is given by Silverstein and Duke^ follow- 
ing AbrikosovS Although the final results will be given 
at T = 0, we work in the space of (fermionic) Matsub- 
ara frequencies for technical simplicity. The rules are the 
following 

(i) In order to get the contribution at the mth order 
in the Kondo interaction Jk, i.e., (^[(Jk)" 1 ] , we draw all 
the topologically distinct diagrams linking m vertices. 

(ii) Each conduction electron propagator (|G12I) is 
represented by a solid line 

fci,iw m e,(7i 



(iii) Each impurity spin pscudofcrmion propagator 
T51 is drawn as a dashed line 
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ECT,CT' 
(20) ((O,{O)=- 



0) 

a' 



(3,0)a(».0>') : 



CO 

a 



CO 

a' 



FIG. 3: Feynman diagram for the conduction electron self- 
energy at the second order in a perturbation theory in Jk. 



Zct,ct' r 
(3 0)b (CO,CO ) : 



l^m* J 71 



(iv) Each vertex 



FIG. 4: Feynman diagrams for the conduction electron self- 
energy at the third order in a perturbation theory in Jk- 



iw m > , 7i 
contributes by a factor 




Jk 



£ 

A=l 



f2,fl 72.71 



(17) 



(v) The independent internal momentums, frequencies 
and spin indices have to be summed over. 

(vi) Each conduction electron loop gives a factor (—A"). 

(vii) Each pseudofcrmion loop contributes by a factor 

(-!)■ 

(viii) Each diagram O [( Jk)™ 1 ] exhibits a combinatorial 
symmetry factor f sym - This factor results from the non- 
compensation between the number of permutations of 
the vertex and the contribution (lJmX) due to Wick's 
theorem^ We show it explicitly in Figs. |3 0] an d HI 
when its value is not equal to 1. 



B. Contributions of the leading order terms to the 
conduction electron self-energy 

In Figs. 130] and[Slwe show the diagrams which will al- 
low us to compute the (3 scaling function perturbativcly. 
As discussed in Section llVl these digrams are sufficient 
to derive the critical exponent of the model. The two 
subscripts in each term of the self-energy expansion indi- 
cate the corresponding powers of Jk and A, respectively 
Each diagram is factorized along two main terms. The 
first one is the spin factor, it results from the summation 
over the internal spin degrees of freedom. In Appendix 
[D]we present the method of calculation of these spin fac- 
tors for all the diagrams we consider. It requires some 
knowledge about the properties (detailed in Appendix 
|A"|) of the SU(A) generators {r} and {<r} associated 
with the impurity and the electronic spins, respectively. 
The second quantity is obtained by summing over the 
internal momentum and frequency degrees of freedom. 
Due to the presence of the imaginary chemical potentials 



{ifj, n }i<n<d {r} , this part has to be computed carefully. 
In particular, the free propagators associated with the 
conduction electron and impurity spins have to be rede- 
fined, which is done in Appendix [U] In Appendix IE! we 

detail the calculation of the third order diagram S(3o) a 
displayed in Fig.0] The other diagrams can be computed 
in the same way. 

The single-impurity Kondo model can be used for the 
study of real systems of dilute (low concentration) impu- 
rities. For that purpose, we adopt the random impurity 
averaging^ the principle of which is the following. A 
number N of impurity spins arc randomly distributed in 
the system. We only keep the linear contribution in the 
impurity density n, = Ni/J\f s . Since the interactions be- 
tween impurities are completely neglected, we recover a 
single-impurity Kondo model. To summarize, each of 

the contributions to the conduction electron self- 

energy corresponding to a diagram of Figs. |3 01 or can 
be written factorized expression as follows 



( _i)«w . { - K y ■ (n^ ■ u v 



x Spin 



■ F< 



<J,l)a 



(18) 



where b p f and I are the numbers of pseudofermion and 
conduction electron loop(s), respectively, Ni is the num- 
ber of impurities in the system, f sym the combinatorial 



symmetry factor of the diagram, Spin S(j^ a the spin 

term corresponding to the contribution of the spin op- 
erators which involves the summation over the internal 
spin degrees of freedom, and F(j^ a the factor resulting 
from the contribution of the fermionic propagators which 
is obtained by performing the summations over the in- 
ternal momentum and Matsubara frequencies. 

We shall exclusively work at T = 0, allowing us to 
calculate the integrals in energy analytically as well as 
the summations over imaginary chemical potentials in- 
cluded in the impurity spin free propagator l|C13|) . We 
follow the conventional cut-off scheme: a cut-off in en- 
ergy is introduced for the conduction electron band, i.e., 
—D < Ek < D, with a density of states p per spin and 
per channel assumed to be constant. 
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At zero-temperature after analytical continuation to the real axis, we obtain the contributions of the first pertur- 
bative order terms to the imaginary part of the electronic self-energy 



Im 



^(2,0) ( W ! W 



(N + 2S- 1)! 2S 
{2S)\(q - 1)\(N - q)\ ' 
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de\- 
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Im 
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(w-ei+iO+)(ei-e a ) 
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hD 
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r}J 
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x Im<J / de 1 de 2 de 3 7 \ ' 7 , - n +u — r ' 6 ( w ~ w 

l- D (w - £1 +iO+)(w - £1 +e 2 - £3 +i0+)(e 2 - £3) 



(20) 



(21) 



where npu is the Fermi-Dirac distribution function. 

In a pcrturbativc RG analysis, we are only interested 
in the energy range < lu <C D. Consequently in the pre- 
vious expressions (fT$)l . lt2U|) and (|2T)l the terms in powers 
of (lu/D) are neglected, and we only keep the factors gen- 



(4,l)a( tn ' lB )= 2 X 



(4,l)b(<M>) = 



Xa.a 
{41)c ((a,(B) = 




crating logarithmic contributions in ln(o;/D), which are 
sufficient to derive the two leading order terms of the 
(3 scaling function as discussed in Section TlV Al When- 
ever necessary, we keep only the principal part of the 
denominator exhibiting a pole. The integrals in energy 
are computed in this limit following the method detailed 
by Ganj2^ we get 



1 



f f + ° 
Im< / de\ — 

I J-d V-£i+iO+) 



-7T , 



"Fd(£2) 



deids2de3 



(£2 - £3) 



nFD(-£l)n F D(£2)"FD(-£3) 



+ "FD(£l)"-FD(-£2)n F D(£3) 



(22) 



Im< / de\de2 , 

[J-d (w-ei+i0+)(ei-e a ) 

7rln(w/D) , 

r+D 

Im > 

_ D (lu - ei +i0+)(w - ex + s 2 - £3 +1O+) 
1 



(23) 



= ir\n.{w/D) 



(24) 



FIG. 5: Feynman diagrams for the conduction electron self- 
energy at the fourth order in a perturbation theory in Jk- 



By incorporating the integrals (1221) , 1|23|) and (|24|l into 
Eqs. 0, 123 and (EH we obtain 
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Im 



(TV + 25-1)! 



(2S)!(g-l)!(iV-g)! 



2S(N - l) - g(iV+ l) 



2S_ 

In 



d m - 1 



{Jk) 2 p- 



5L 



N 



(j K p) ln(u / D) + ^(j K p) 2 K ■ ln(u / D)[ . (25) 



IV. PERTURB ATIVE RENORMALIZATION 
GROUP APPROACH 

In this Section we use a perturbative RG method, the 
Anderson's "poor man's scaling" approachi2£i2£ The in- 
termediate coupling (IC) fixed point associated with the 
model is identified, for various limits characterized by 
the shape of the {r} representation. We start by deriv- 
ing the (3 scaling function, as a function of ( Jkp) ■ A fixed 
point corresponds to a zero of the /3 function. A simple 
analysis of the scaling function properties around a fixed 
point allows also to determine its stability**! In the case 
of a stable IC fixed point, a rich physics characterized 
by a non-Fermi liquid behavior of the system (namely, 
anomalous temperature dependence of the specific heat, 
magnetic susceptibility ...) can be expected. 



A. Derivation of the scaling function from the 
conduction electron self-energy 



In the previous Section IIII Bl we have calculated the 
first order terms of the conduction electron self-energy in 
a perturbative approach. This yields the leading terms 
of the (3 scaling function 

d(J K p) 



P(Jkp) 



d(\nD) 



(26) 



The connexion between both results from RG arguments, 
more precisely the Callan-Symanzik cquatioi 



,68.69 



as ex- 



plained below. 

Without mentioning the spin dependence, we rename 
the imaginary part of the conduction electron self-energy 



Im 



J {r> 



, Eq. 



J {iT 



to have shorter notation 



in what follows. Following Gan et al.M*z± S^ r j can be 
expressed in the most general way up to the fourth order 
in perturbation theory as 



J {r}t 



D, J K p) 



(27) 



Po{Jkp) 2 + PiHuj/D)(j KP y 



P 2 ■ K -\n(cj/ D)(J KP ) 
P 4 \n 2 (uj/D)(j K p) 4 + 



p 3 -K-(j KP y 



(Jkp)' 



where a, Po, ... , Pa are coefficients independent of lu, 
-D, JkP- P^{r} i s the dimensionless scattering rate (i.e., 
inverse of the quasiparticle lifetime). The key point of 
the RG theory is that the underlying physics must be 
independent of the renormalization scheme, namely the 
cut-off D. Therefore the scattering rate is invariant un- 
der RG transformation. The Callan-Symanzik equation 
resulting from this property is then 



d 



d(hxD) 



+ P(Jk P ) 



d(J K p) 



pE {r} (w,D, J K p) 



= . (28) 

The /3 scaling function is taken in its most general ex- 
pansion up to the third order in powers of JkP 

(3{JkP) = ni{j K p) 2 + k 2 ■ K ■ (J K p) 3 

+ n 3 (j K p) 3 + o[{Jkp) 4 ] ■ (29) 

Eqs. l(2T|) and (J2U incorporated into Eq. l(2*5|l give 

[2 Kl P - P x ] (J K p) 3 + [2 K2 P - P 2 ]K ■ (J K p) 4 

+ 2 K3 P (J K p) 4 + [S^Pj - 2P 4 ] ln(w/D) • ( Jkp) 4 



O 



(Jkp) 1 



= 



(30) 



Eq. (|30|l generates logarithmic series valid for any cou- 
pling constant Jk- Therefore the coefficients in front of 
the terms of successive orders in Jkp are identically 0. 
The result is 

(3D 



K 2 
Pi 



2P 

2P ' 
, 

2 



(32) 
(33) 
(34) 



The coefficients Po, Pi and P 2 are obtained from the 
imaginary part of the conduction electron self-energy, Eq. 



P> 
Pi 
Pi 



1 

N ' 
-1 , 
1 

2 ' 



(35) 
(36) 
(37) 
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We point out here that the coefficients given by Eqs. (|35|l . 
(|36[) and (|37|> are exact at this stage for any value N, 2S 
and q, without any restriction. As far as the consistency 
equation 1|34|) is concerned, it has been checked by Gan2£ 
for an impurity spin in the fundamental representation 
of the SU(2) group. At each order of perturbation our 
study leads to the same integrals in energy displayed in 
Eqs. (23) an d 123) • The only difference between 

Gan's work and our generalized approach is contained in 
spin factors which depend on the considered group and 
representation describing the impurity. Since we have 
calculated these spin factors in an exact way, we conclude 
that the previous consistency equation is automatically 
satisfied for any group and any impurity representation 
we are interested in: the explicit calculation of P4 for the 
generalized {r} representation is not required. 



B. Analysis of the intermediate coupling fixed 
point 



the expansion of the (3 function, O (Jkp) 4 , contains a 
coefficient proportional to the eigenvalue of the quadratic 
Casimir operator IjAllf) related to the impurity spin rep- 
resentation. In order that the latter contribution which 
has not been taken into account here remains negligible 
compared to the second and third order contributions, 
we must also suppose that N, 2S and q are small com- 
pared to K. To summarize, the domain of validity of the 
scaling function perturbativc expansion is given by 



A > 1, N < K, 2S <C K, g « A 



(40) 



Consequently the IC fixed point describes an over- 
screened impurity, as discussed originally by Nozicres and 
Blandini The single-channel case turns out to be out of 
reach of the present study. 

In addition to the value of J*, another key quantity 
which can be drawn from the (3 function is its slope A at 
the IC fixed point, namely 



Following the results of previous Section TlV Al the (3 
scaling function (|26|) obtained in a perturbative theory 
calculation in Jk up to the third order is given by 



P(Jkp) = 



N, .a 



NK 3 
-^-( j kP) 



O 



CM 4 ' 



(38) 



One of the main results of this work is that Eq. i|38|) is 
valid for any value of the parameters 25 and q character- 
izing the L-shaped representation of the impurity spin, 
provided that the perturbative expansion in Jk is valid 
as discussed in the next paragraph. 

The critical value J* taken by the Kondo coupling at 
the IC fixed point is obtained from the zero of the (3 
function, i.e., /3(J*p) = 0. We get 



J P 



2 

A 



O 



1 



(39) 



The latter expression for J* allows one to determine the 
domain of validity of the perturbative expansion. The 
RG perturbativc approach we use is valid in the weak 
coupling regime only. The IC fixed point has to be in 
the vicinity of the free electron fixed point of the sys- 
tem. Consequently, the perturbative expansion of the (3 
function is trustworthy exclusively when J* is close to 
0, which means by using Eq. I|39|l K ^> 1. Since each 
conduction electron loop which gives a factor of (—A') 
contains two interacting vertices at least, and each term 
J*p contributes by a factor of 1/K, we conclude that 
the diagrams displayed in Figs. |3 El and include all the 
contributions to the electronic self-energy up to the order 
O [1/ A 3 ] at the IC fixed point. This argument allows one 
to check the self-consistency of the calculation. Another 
restriction is due to the domain of validity of the pertur- 
bative expansion. As previously noticed by Abrikosov- 
MigdalZS and Affleck-LudwigfS the fourth order term in 



dd 



d(J K p) 



Jk=J* 



N 
A 



O 



1 



(41) 



In the limit given by Eq. H4U|) that we consider, A is 
automatically positive therefore the IC fixed point is 
stabled The slope A dctc rmines the critical exponent 
of the model. It is a strong mark of the universality con- 
cept characteristic of the Kondo effect in the low-energy 
regimeJi*££ This regime describes the system at temper- 
atures that are low compared to the bandwidth D. In 
that case the physics of the system is governed by a sta- 
ble IC fixed point, and the temperature dependence of 
the various physical quantities at low T obeys a power 
law in T/Tk which is entirely controlled by A (7k being 
the Kondo temperature). 

The Kondo temperature is the only relevant energy 
scale in the problem, it is generated dynamically by the 
interactions with the impurity spin. It plays an essential 
role in the Kondo effect: 7k characterizes a crossover 
between the weak coupling regime and the low-energy 
regime. This temperature scale is invariant under the 
scaling process. One can easily derive it by inverting and 
integrating the (3 function expression^ By using Eqs. 
(2I>J) and (211) we get 



k B T K ~ D(J K p) K ' N exp 



N(J K p) 



(42) 



As mentioned above, the results we have obtained arc 
valid in the overscreened regime for any value of the pa- 
rameters N, 2S and q, provided that Eq. H4()(l is satisfied. 
The impurity spin representations that are adressed by 
the present study are displayed in Fig. HO We conclude 
that all the Kondo models defined by these different im- 
purity spin representations belong to the same universal- 
ity class in the overscreened regime, which is character- 
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(a) 



□ 



(b) 



(c) 



(d) 



FIG. 6: Irreducible SU(iV) representations of the impurity 
spin considered in the present work: (a) L-shaped ; (b) fun- 
damental, 2S = q = 1 ; (c) symmetric, q = 1 ; (d) antisym- 
metric, 2S — 1. 



ized by the following low-temperature behaviors 2 ^ 
— resistivity due to Kondo scattering : 



S Pe (T) = p e (T = 0) - a 



T 



impurity specific heat : Ci mp (T) 



magnetic susceptibility : Ximp{T) 



T 



2A 



1 / T 



TVTi 



2A 



We can notice that the system described by an IC fixed 
point exhibits non-Fermi liquid properties. As in the 
original single-channel Kondo model characterized by an 
exactly screened spin 1/2 impurity, the resistivity in- 
creases when the temperature is lowered. Denoting by 
S e ff the effective spin of the impurity screened by the 
conduction electrons and ( ) the thermal average, we 



have^ ((S eff ) ) ~ T Ximp (T) ~ T 2A as T 

which indicates that the IC fixed point is associated 
with a ground state corresponding to a spin singlet. In 
the overscreened regime the impurity spin is completely 
quenched by conduction electrons. 

The value of the universal critical exponent A we have 
obtained in Eq. I|41|) is in agreement with the result 
of Gan et at who considered a spin 1/2 in the SU(2) 
g r0U p ( 24 i 25 r-Q^ cage Q £ i ar g er S p ms > 1/2) in the 

SU(2) group has also been studied^S^i leading to sim- 
ilar results. Eq. (|41|) allows one to recover the over- 
screened limit of the exact solution of the multichannel 
Kondo model, obtained by using the Bethe ansatz for 
the symmetric impurity spin representations of SU(A f ) ) 9 . 
and boundary conformal field theory in the case of an ar- 
bitrary impurity spin representation of the SU(2)5£ and 
SU(iV)S groups. It also corroborates the large- N stud- 
ies of the SU(./V)xSU(.K') Kondo model involving either 
the fundamental^ or the antisymmetric^ impurity spin 
representations, although the K ~ N limit considered in 
these works is different from ours. 



V. SUMMARY AND CONCLUSIONS 

In this paper, we have studied a multichannel SU(iV) 
Kondo model, involving a mixed symmetry impurity spin 



representation. Namely the generalized impurity spin is 
realized by a combination of 2S bosonic and q fermionic 
degrees of freedom, in the presence of K conduction elec- 
tron channels. The study has been carried out by making 
use of a perturbative renormalization group technique, 
the Anderson "poor man's scaling" approach. The main 
difficulty for the description of spin systems is that spin 
operators obey neither Fermi nor Bose statistics, leading 
to the absence of Wick's decomposition theorem in this 
context. Consequently conventional perturbative meth- 
ods are prohibited. We have used the Abrikosov pseud- 
ofermion formulation to describe the impurity spin states, 
and tackled the problem by making use of a method initi- 
ated by Popov and Fedotov. Its basic principle consists in 
introducing a set of auxiliary chemical potentials which 
eliminates the unphysical Hilbcrt sub-space keeping the 
physical spin states only, and therefore enables a proper 
handling of the fermionic constraint. 

The partition function has been obtained within a path 
integral approach. The general expression of the conduc- 
tion electron self-energy imaginary part has been derived 
perturbatively in powers of the Kondo coupling constant, 
i.e., by considering the weak coupling limit. At each or- 
der of a perturbative calculation the derivation involves 
two steps, namely the calculation of the spin factor and 
the computation of the fermionic propagator contribu- 
tion which includes the Popov-Fedotov auxiliary chem- 
ical potentials. The imaginary part of the conduction 
electron self-energy is proportional to the electronic scat- 
tering rate, which is a physical quantity and therefore 
is invariant under renormalization group transformation. 
The Callan-Symanzik equation can be applied, and yields 
the derivation of the two leading terms of the (3 scaling 
function, i.e., up to the order C[l/i^ 2 ]. We have estab- 
lished that this truncated expansion of the (3 function is 
valid in the overscreened (Nozieres - Blandin) regime, for 
any value of N, 2S and q providing that these parameters 
remain small compared with K. The intermediate cou- 
pling (IC) fixed point describing the low-energy physics of 
the model has been identified in the overscreened regime. 
The critical exponent A, corresponding to the slope of the 
p function towards JkP at the IC fixed point, is given by 



A 



N 
K 



O 



1 



A governs the temperature dependence of the various 
physical quantities at low T which are of non- Fermi liquid 
type, as discussed in Section IIVI The IC fixed point 
describes the system in its ground state which is a spin 
singlet: in the overscreened regime the impurity spin is 
completely quenched by the conduction electrons located 
on its site. 

Several types of impurity spin representations have 
been considered in previous works and can be viewed 
as special limits of the mixed symmetry formulation we 
have used. Hence direct comparisons with our study have 
been made possible. Our results are in agreement with 
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those obtained by other methods, like Bcthc ansatz and 
boundary conformal field theory. We have established 
that the Kondo models defined successively by funda- 
mental, symmetric, antisymmetric and mixed symmetry 
impurity spin representations yield exactly the same low- 
energy physics in the overscreened regime. To summa- 
rize the present work gives a unified view of the over- 
screened Kondo regime, for a large set of representations 
of the impurity spin. It is also possible to consider ar- 
bitrary impurity spin representations, beyond the cases 
mentioned here. For example one can think about repre- 
sentations of the SU(/V) group with rectangular-shaped 
Young tableaux, as discussed in Ref. 0. The method 
we present can be generalized to study the case of an ar- 
bitrary SU(1V) representation. No additional conceptual 
difficulty has to be expected. The spin factors associated 
with such a representation, which arc present in the con- 
duction electron self-energy, have to be computed follow- 
ing the derivation displayed in Appendix [D] The traces 
over r matrices we have to consider involve products of 
two or three generators in the impurity spin representa- 
tion, and therefore the spin factors can be derived ex- 
actly for any irreducible representation of SU(N)m££l£^ 
However the computation of these spin factors as a func- 
tion of the Young tableau parameters describing a com- 
pletely arbitrary impurity spin representation certainly 
requires extremely lengthy algebra, as the complexity of 
Eqs. dD9j, G5B and lfDT3|i . which are restricted to the 
L-shaped case we consider, suggests it easily. 

Finally we mention the case of the single-channel 
SU(-/V) Kondo model involving a mixed symmetry im- 
purity spin, which we considered in a previous work4^£ 
In such a case the L-shaped impurity exhibits a rather 
surprising duality: it is exactly screened if one considers 
the fermionic component, but underscreened towards the 
bosonic one. An interesting problem remains. Namely, 
the understanding of the low-energy physics in the regime 
characterized by an instability of the strong coupling 
fixed point turns out to be out of reach of the present 
work. This issue deserves further investigations. It might 
have important future applications for the Kondo lat- 
tice problem^ with potential implications for the de- 
scription of heavy-fermion compounds. In order to gain 
some insight it would be interesting to carry out a non- 
perturbative study. Unfortunately it is quite a difficult 
task. Indeed well-known methods usually efficient to 
study the Kondo problem, as numerical renormalization 
group or Bethe ansatz techniques, have severe limitations 
in the context of the mixed symmetry representation. For 
example, the simplest impurity that displays an IC fixed 
point in the single-channel case implies a large Hilbcrt 
space which restricts the efficiency of a numerical renor- 
malization group study. The difficulties encountered by 
the Bcthc ansatz are of different nature: due to the mixed 
symmetry properties of the impurity spin representation, 
the model with Kondo coupling only is not integrable^4 
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APPENDIX A: PROPERTIES OF THE SU(iV) 
GENERATORS 



In this Appendix we give some properties of the gener- 
ators {cr} and {r}, in the case of the two types of irre- 
ducible representations we have to consider in the SU(iV) 
group. The first one concerns the conduction electron, it 
is the fundamental (single-box) representation {cr} dis- 
played in Fig. |5J Its dimension is 



d {a} = N, 



(Al) 



and the eigenvalue of the quadratic Casimir operator cor- 
responding to {cr} is 



c 2 (M) = 



N 2 



1 



2N ' 



(A2) 



which in the case of a spin S in SU(2) takes the usual 
value S(S +1). The associated traceless hermitian ma- 
trices {&" A } 1<A<N 2_ 1 obey the following properties^ 

— multiplication law: 

~a~b 1 <f f 
cr cr — diRif^i 



N'-l 



N 2 -l 

(d ABC 
C=l 

— normalization relation: 

~Z$AB , 



if 



ABC) a C 



cr A cr B ] 



Tr[ 



completeness relation: 



(A3) 



(A4) 



N -1 

E„A A 



1 



8 a ,s8p n - jy6 a> p5y t s ) , (A5) 



where A, B and C label the SU(iV) generators matri- 
ces (1 < A, B, C < N 2 - 1), J w is the (N x N) iden- 
tity matrix, a, (3, 7 and 6 label electronic spin colors 
(1 < a,/3,7,(5 < N), and d ABC and f ABC are the struc- 
ture constants related to the d and / Gell-Mann tensors, 
respectively^ The structure constants are characteristic 
of the underlying Lie algebra. These tensors follow the 
rules 

- d is totally symmetric: d ABC = d BAC , 

- / is totally antisymmetric: f ABC = -f BAC , 

- d and / are traceless tensors: 

Y,d ABB 



f ABB = 



B 
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In SU(2), d = and f = e, the totally antisymmetric 
(Levi-Civitta) tensor of rank 3. In further calculations, 
we use the following idcntitiesB*^ 

J2 d ACB d GBC = ——^—-5 ag , (A6) 

B,C 



cGBC 



J2f ACB f 

B,C 

J2d ACB f GBC = 

B,C 



N 
NS AG 



(A7) 
(A8) 



The second representation of interest deals with the 
impurity spin. It is a L-shapcd representation denoted 
by {r}, and characterized by 25* boxes in the first row 
(horizontal line) and q boxes in the first column (vertical 
line), as displayed in Fig.^ Its dimension can be derived 
by the use of Robinson's formula^ 



25 



Z{r} 



n 2S r>?-i 

~ L, Af+2S-l°jV-l J 



(A9) 



25* + q 

and the eigenvalue of the quadratic Casimir operator as- 
sociated with {r} is given bjiS 



c 2 ({r}) 



1 f {2S + q - 1)[N 2 - {2S + q - 1)] 



N 



N 



(A10) 



where nij is the number of boxes in the j-th row (1 < 
j < N) of the Young tableau describing {r}. Finally 



c 2 ({r» 

25 + q - 1 



2S-q + N- 2S + « r ) 



We give now useful identities satisfied by the (rf{r> xc ^{r}) 
components of the {r^ 1 } matrices. The normalization 



relation 



,80 



^{r}^({r}) S AB 
~~ d W} C 2 ({*}) 2 



(A12) 



The generators in the {T} representation obey also some 
particular multiplication laws for the contracted indices 
close to each other— 

^f B f A f B = |c 2 ({r}) - C2 (^ dj » |f^ ; ( A13) 



J2t c t B t B t A = C 2 ({T})t c f A , 



(A14) 



where C2({Adj}) is the eigenvalue of the quadratic 
Casimir operator associated with the adjoint representa- 
tion of SU(AT). Due to the choice of normalization given 
by Eq. (JUJl we have: C 2 ({Adj}) = N. 



APPENDIX B: 



POPOV AND FEDOTOV 
METHOD 



In this Appendix we briefly explain a method pio- 
neered by Popov and Fedotov^ in order to handle the 
constraint JHJ) fixing the number of pseudofcrmions on 
the impurity site. Originally formulated for spins 1/2 
and 1 of the SU(2) group, this powerful method has 
been developed for arbitrary spin in SU(2)2£ and gener- 
alized to SU(7V)— It has already been used to study the 
Heisenberg model^i the Kondo latticed and disordered 
systems—^ 

The principle is to express the partition function of 
a spin system in terms of partition functions of a corre- 
sponding fermionic system, pseudofermions being used to 
describe the spin states. The problem is to incorporate 
the fermionic constraint properly, namely to eliminate 
the unphysical Hilbcrt sub-space without removing the 
true spin states. This can be done by noticing that the 
Kondo hamiltonian (jSJ without imposing any constraint 
is disconnected in the peudofermion charge sectors. Con- 
sequently the trace over the full Fock space (i.e., without 
constraint) can be rewritten by applying the Pauli prin- 
ciple as 



Tr Fock exp 



f3(H K - fihf) 



Tr n/ exp - /3(H K 



n,=0 



(m f ) 



(Bl) 



where Tr n/ is the local subtrace over the states with a 
fixed number nf of pscudofcrmion(s) on the impurity site. 
The physical partition function Z <|1U|) . which includes 
the fermionic constraint, is together with Eq. (|B 1|) 

Z = 5{rif - l) • Tr Fock exp - /3(7Y K - Mn/) . (B2) 
By using the identity 

8{n f - l) = — — V e iT(2«-l)/d{r} e iir(X-2n)n//d { r} ) 
n— 1 



incorporated into Eq. (|B2|) we obtain 



Z = j Vfi P{fx] Tr Fock cxp - p(Hn - fihf) , (B3) 
with 



»{r} 



1{T] 



-(2n - 1) 



X *l" + ^ 2n - 1 )) ' (B4) 



which is exactly Eq. Ijllfl. 
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APPENDIX C: DERIVATION OF THE FREE 
PROPAGATORS 

In order to compute the conduction electron self- 
energy at the first orders in a pcrturbative theory in Jk, 
the first step consists in deriving the free propagators 
associated with the conduction electrons and the pscud- 
ofermions. It is performed in this Appendix. Due to 
the presence of the imaginary chemical potentials l|llfl in 



the partition function Z i|16|) . we should determine the 
free propagators carefully. For that purpose we use the 
technique of the generating function^l^ 

Let us denote by J*, J and 77*, 77 the external sources 
coupled to the conduction electron c, c* and impurity 
fermion /, /* variables, respectively. These external 
sources are Grassmann numbers. The generating func- 
tion is 



Z G [3,r,r,,ri*] 



(CI) 



f ( 1 d{r} 

/ T>c"DcDf*'Df\ V" exp [i/?/j„] • exp - V" V" {c£ a x (iu m )(-iu) m + £fe)cfe, CT ,\(iui m ) 



_ Jfe, CT ,x( ia; m)Cfe, < T,A(iw TO ) - Cfe . A (iw TO )Jfe i0 . j A(iw m )} 

E E Uj( iuJ m )(-uj m + i/i„)/ 7 (icj m ) 



-?7*(iw m )/ 7 (iw m ) - /*(iw m )?7 7 (iw m )} 



If each source is put to one can see immediately that Z G is identical to the free partition function Zq, describing 
the system in the absence of the Kondo interaction. The free propagators are given by the thermal averages 



(Ck ,<r ,A ( iw "io )4 ,<r , A ( iu "o )) 

(/ 7o (iw mo )/* o (iw mo )) 



1 &Z<>[J,J\r,,T,*] 
Z dJ ko ,a ,\ (i^ro )9Jfe OlO - ,A ( iw ™o ) 

1 d 2 z G [j,r, v , v *} - 

Z Q d?] 10 (iuj mo )dr)* (iuj mo ) J J = J* = 

T) — TJ* — 



J — J* — 
r) — rf — 



We use then the following identity for Gaussian integrals over Grassmann fields 
/ Vc*VcDf*Vfexp[-J2 E {3Ua(^™)( 

-iw m + £fc)c fei ^ A (iWm) 

iuJ m k,cr,X 

- EE {/ 7 (i^m)(-iw m + i//„)/ 7 (io) m ) 



exp 



-77*(iw m )/ 7 (iw m ) - f*(iLL> m )rj 7 (iuj m )} 

E E k^w — 

-1 



h 7 (iw m ) 



x / Vc*VcDf*Vfcxp 



\LO m — \\X n 

E E {Cfc : <r,A( iw m)(- iw m + £ fc)cfc,<T,A(iw m )} 



icj m k.a.X 



^2^2 {f^ XU) rn){-iU m + i/X n )/ 7 (i^m)} 



(C2) 
(C3) 



(C4) 
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By using Eq. (|C4|) one can easily obtain the derivatives of the generating function Z G (IClll with respect to the sources, 
allowing us to get the thermal averages 



1 



d 2 Z G [3,J*,ri,ri*] 



Z dJ fco!CTO ,A (i w mo)dJfc 0:CTo Ao (iw mo ) 

1 d 2 Z G [J,J*, V ,r,*} 



1 



J = J* = 

Tj — TJ* — 



Z d-q 10 (iw mo )dri* (iw mo ) 



J = J* = o 
rj — 77* — 



d 



{r} V 



mo ^fco 
1 



1^ Ee^^(i^ 



(C5) 



(C6) 



where Zq is the free fermionic part of the partition func- 
tion, associated with the imaginary chemical potential 
(i/^n) 

Z&M (C7) 

= fvf*Vfcxp[- {/*(^m)(-iw m + i/Xn) 

iu m 7 

x / 7 (iw m )} 

In Eq. l|C7p we recover the partition function of free 
fermionsSi of "energy" (i/U n ) 

Z$(Vn) = f( [l + e-V»~] = [l + e-^] d{r > . (C8) 

7=1 

Due to the particular form of the imaginary chemical 
potentials, one can derive some identities with the help 
of simple trigonometric properties, which will be useful 
to perform the calculations involving the pseudofermion 
propagators (see Appendix lE|l 



^e i ^»[l + e- i ^]^ r > = [d {r} ] : 

71=1 



(C9) 
(CIO) 



^e^^l + e-^]^" 1 =d {r} (d {T} -l) . (Cll) 



Using Eqs. (|C2|) and (|C5|) we get for the free conduc- 
tion electron propagator 



(Cfc , CTo ,Ao(iw mo )4 oCTOi A (iWmo))o 



1 



(C12) 



lW mo — Sk 

The expressions l(C6)l. (|U8l and (EU yield the 

free propagator associated with the impurity spin pseud- 
ofermion 



(f~f„(ium )f* (iuJ mo ))o (C13) 



[d 



mi „=i 



APPENDIX D: COMPUTATION OF THE SPIN 
FACTORS PRESENT IN THE CONDUCTION 
ELECTRON SELF-ENERGY 

In this Appendix, we present the method of calcula- 
tion of the spin factors necessary to obtain the electronic 
self-energy perturbativcly. The properties of the SU(A^) 
generators given in Appendix ^ will be used intensively. 
We will detail the contribution of each diagram of Figs. 
[3J 0] and [3] A general remark one can make is that the 
impurity spin propagators form a closed loop in each di- 
agram. Therefore the {t* A } matrices will always appear 
as a trace over products of generators in the {T} repre- 
sentation. 



1. Spin multiplicative factor in the second order 
term of the perturbation calculation of the 
conduction electron self-energy 

After summation over the internal spin degrees of free- 
dom, the spin factor contribution of the second order 
term to the electronic self-energy is given by 



Spin 



J (2,0) 



jv- 



£ Tr[f^f B ] • [a B a\, a . (Dl) 



A.S = 1 



The SU(A^) generators in the {r} representation obey 
the following normalization relation^ 



Tr[r A r B ] 



^ C2 ( {r} l T r[^l 



d {a} C 2 {{(r}) 



(D2) 



which gives by making use of Eqs. (|A1|) . (|A2|> . HA4jl. (|A9|) 
and !jATU|> 



Tr[f A f s ] 

(N + 2S-1)\ 



(D3) 



25 



{2S)\{q - 1)\(N - q)\ 2N(N 2 - 1) 

x [2S(N - 1) - q(N + 1) + (A^ 2 + l)]6 AB , 

where 2S and q are the parameters of the Young tableau 
associated with the L-shaped representation, as shown in 

Fig.m 
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Using Eqs. I|D3|) and (|A5|) wc obtain finally the spin 
multiplicative factor of the second order term 



Spin 



J (2,0) 

(N + 2S 



1)! 



25 



(2S)\{q - l)\(N - q)\ AN 2 

x \2S(N - 1) - q(N + 1) + (iV 2 



(D4) 



In the previous expression the summation is performed 
over all the permutations of the generator indices. 



STr 



ACB 

{r} 



can be calculated for any representation in 



SU(iV) (N > 3) by following the method developed by 
Okubo^ that we summarize here. The Young tableau 
in the standard arrangement of the L-shaped represen- 
tation is specified by TV integers (/i)i<i<jv, fi being the 
number of boxe(s) in the ith row 



2. Spin multiplicative factors in the third order 
term of the perturbation calculation of the 
conduction electron self-energy 



/i=25 

fi = 1 , 2 < i < q 

fi = , q + 1 < i < N. 



After summation over the internal spin degrees of free- 
dom, the first spin contribution of the third order term 

is given by the diagram corresponding to S/^\ a 



Spin[E^ ff 0)a 

N 2 -l 

= Yl Tr[f^r c f B ] • [* B * c & A ] a , c 



(D5) 



A,B,C=1 

The calculation of Tr[f A f c f B ] for the L-shapcd rep- 
resentation in the SU(iV) group is rather technical. For- 
tunately this kind of algebraic object has already been 
intensively studied in the field of particle physics. We can 
express it with the aid of a symmetrized trace, which de- 
fines a symmetric invariant tensor associated with 



Tr[t A t c fB] = STi acb + i rf { r}C 2 ({r}) / 



ACB 



2 A 2 - 1 

where STr^j 5 is the third order symmetrized trace 



, (D6) 



STr AiMA 3 



3! 



Tr[f A "< i >r A *( 2 ir A »(3>] 



Wc define the indices (<7j)i<j<N by 



then 



with 



N + l 



1 N 
A^ /? 



N + l 2S + q-l 
~ A 



0" - /,-), 



STr 



ACB 

{r} 



d ACB 



K{{T}) 



K({T}) 



N 



N 



{N 2 -l)(N 2 -i) 



(D7) 



(D8) 



After lengthy algebra we get using Eqs. (|A9() . (| A10|> . 
(lD6ll . lHT7tl and 



Tv[f A t c f B ] = 



(AT + 25- 1)! 



25 



(2S)\{q - 1)\(N - q)\ 2{N 2 - 1) 
f 1 



(N 2 - 4) 



(2sy 



2 

A" 3 



A 



— + 3 + A 



+25 



+ 



25 



2 

A 

1 " 

N 



4 

7N A 3 
— + — 

q N 



7N 3N 2 



2 2 

d ACB 

N- 



+ 2S-q 



4 

N' N 



A 



7N 3N 2 
~~2 2~ 



cACB 



A 



(D9) 
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We have checked the validity of the previous expres- 
sion (|D9(1 by comparing it with the results of Okubo 
and Patera*^ who considered special L-shaped represen- 
tations, i.e., with given value of 25 and q for general 
SXJ(N) groups. We get a perfect agreement. 

For the electronic part we can use 



\a B a c a A ] , 

L J a c 



1 



2A 



(D10) 



G=l 



By combining Eqs. I|A5[) . l|A6jl . (|A7|) . 

previous results (|D9|> . ijDIOf) . we obtain 



Spin 



(3,0)a 

(A + 25-1)! 



81 and the 



(Dll) 



25 



(25)!(g-l)!(JV-g)! 87V 3 

x [ (25) 2 (2-3A + A 2 ) +<7 2 (2 + 3A + A 2 ) 

+2S ■ q(4 - N 2 ) + 25( - 4 + 3A - AN 2 + 2A 3 ) 
+q( - 4 - 3N - AN 2 - 2N 3 ) 

+ (2 + 4A 2 + A 4 ) 



The second spin contribution of the third order term 



corresponds to the diagram S 



(3,0)b 



Spin [E^' 



0)1., 



(D12) 



= £ Tr[f A f c f S ] • [& c & B a A ] a , c 



A,B,C=\ 



which can be derived in a similar way. The result is 



Spin 



(3,0)b 

(A + 25-1)! 25 



(D13) 



(2S)\(q - 1)\(N - q)\ SN 3 

(25) 2 (2 - 3A + N 2 ) +q 2 (2 + 3N + N 2 ) 
+25 • q(4 - N 2 ) + 25( - 4 + 3A - 3N 2 + N 3 ) 
+q( - 4 - 3N - 3N 2 - N 3 ) + (2 + 37V 2 ) 



3. Spin multiplicative factors in the fourth order 
term of the perturbation calculation of the 
conduction electron self-energy 

The first spin contribution resulting from the fourth or- 
der term is given by the diagram corresponding to S^"^ 



Spin [E^" 



(D14) 



N — 1 



J2 Tr[r D f c T B T A ]-Tr[a B v c ] 



A,B.C.D=1 



X \& D * A ] 



With l|AT|) and (|A14|I we get 

J2TT[t D t c t B t A ] -Tr[& B & c ] 



(D15) 



Bfi 



(A + 25-1)! 25(25 + ^-1) 



(25)!(g-l)!(A-g)! 8A 2 (A 2 -l) 
x [25(iV - 1) - q(N + 1) + (A 2 + l) 

Using Eq. (El we have 

Spin [^ (4 ,l)a 

(A + 25-1)! 25(25 + q-l) 



Sad ■ 



(D16) 



(25)!(g- l)I(JV-g)! 16A 3 



25 (A - 1) - q(N + 1) + (A 2 + 1) 



The second spin contribution resulting from the fourth 



order diagram ^ b is 



N z -1 



Spin 



E<7.(7 
(4,l)b 



E Tr[f A f B ] -Tr[f c f c ] 

A,B,C,D=1 

xTr[^].[^] CTV .(D17) 



By using Eqs. HA4f) . (|A5|) and (|A12|> we get 

Spm^L (41)b 

/ (A + 25-1)! 25 
: 1 (2S)l(q - 1)!(A — q)\ ' ~N 

25(A - 1) - q(N + 1) + (A 2 + 1) 
1 



(D18) 



16A(A 2 -l) a - a 
The last spin multiplicative factor showing up in the 
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fourth order diagram S^J\ C is given by 



(4, 
N 2 -l 



Spin 



2j (4,l)c 



Tr[t A t c t D r B ] (D19) 

A,B,C,D=1 



The expressions (IA4|) and l|A13() lead to 

J2TT[t A t c t D t B ] -Tr[a B a c ] 



(D20) 



B,C 



(N + 25-1)! 



25 



(25)!(g-l)!(JV- 9 )! 4A^(iV 2 - l) 
25 + ^ 1 -{25(7V-l)-g(7V + l) 



2N 



N 

T 



+ (7V 2 + 1)} 

x \2S(N -I) - q(N +1) + (N 2 + 1] 
Using Eq. we obtain 



Sad 



Spin[E^ 1)c j (D21) 

(JV + 2S- 1)! 25 _ 

(2S)\(q- l)\(N-q)\ ' 167V3 ' 

x | [25(7V - 1) - q(N + 1) + (W 2 + 1)] 2 

x (25 + g-l) 

25(iV - 1) - q(N + 1) + (N 2 + 1) iV 2 | . 



APPENDIX E: EXPLICIT CALCULATION OF A 
THIRD ORDER PERTURB ATIVE DIAGRAM 

In this Appendix we show how to compute the contri- 
bution due to the diagram associated with Spj^)\ a in Fig. 



0] All the other diagrams can be treated in the same way. 
The spin multiplicative factor resulting from the {T A } 
and \er A } m atrices has been calculated in the previous 
Section TP 21 it is given by Eq. (|D11() . Here we just show 
how its starting point (|D5(I can be derived. Then we will 
focus on the frequency and momentum part. In Fig. [7| 

the diagram associated with S/^q^ is represented, by 
including all the internal degrees of freedom. 

Let us start with the spin contribution. According to 
the rules defined in Section llll Al we obtain by expliciting 
the diagram shown in Fig. [7] 
Spin 



N 



i\r-i 



E E (E<. r 

CTi,(T2 = l 71,72,73 = 1 A=l 

N 2 -l 

■ E 



A 

71,73 



ct',CT2 72,71 



B = \ 
N 2 ~l 

E 

c=\ 



c v c 

a 2,01 73,72 



[^ B a c a A ] a ,^Tr[t A t c t B ] 



(El) 



A,B.C 



which is identical to Eq. i|D5l 



We consider now the contribution including the fre- 



quencies and the momentum, denoted by Fi 



(3,0)a 



Al- 



though our analysis is restricted to the zcro-tcmpcraturc 
limit, the calculations are performed in the space of 
(fermionic) Matsubara frequencies for reasons of simplic- 
ity. We shall use the properties of the free propagators 
presented in Appendix [CJ The rules defined in Section 
IIII Al combined with expressions l|C12ll and (|C13(I give 




After summation over the imaginary frequencies we get 
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^(3,0)a(H,e>iCV') = 



"2 



\ ivy 



FIG. 7: Feynman diagram of the third order in perturbation for the conduction electron self-energy (according to the notations 
defined in Section fill At . 



(3,0)a 



yVs/ fex,fe 2 \Rr}J m / 

(— ^-5 E^Il+e* 1 ]^ 1 J ' (^-2 E c ' 
\[ d {r}] « 2 / \[ d {r}] „ 3 



» 1 + e 



-i/W,l d {r} 



x [n FD (e fel ) - n FD (i//„ 3 )] ■ [n FD (e fc2 ) - n FD (ifx n2 )] 

n FD (i^ m e - e kl + i/in 3 ) ~ ^Fp(iMni) 

(iw TO e — e fel — i/z ni + i/i„ 3 )(efc 1 — £fc 2 + ifi n . 2 — i^„ 3 ) 

riFp(iaw - gfc 2 + 'Mn 2 ) - ^Fp(iMni) 

(iw m e - £ fc2 - i/i ni + i/i n2 )(-e fel + £ fc2 - i/i„ 2 + i/i„ 3 ) 



(E3) 



where n FF j is the Fermi-Dirac distribution function. We see in Eq. 1|E3(I that the poles of the free propagators are not 
on the real axis, as usual. Resulting from the handling of the fermionic constraint 1)1 1)1 . the poles are shifted by the 
imaginary chemical potentials i/i„. By using the following contraction identity 



["-FD^fcJ - 7l FD (1/^13 )] ' [«FD(iw m = - E kl + i/V, ) - JlFD^nj] 

-n FD (-£ fcl )7i FD (-i / u ni )n F D(iMn 3 ) - n FD (£ fcl )n F D(iAini)«-FD(-iMn 3 ) 



(E4) 



we have 



(3,0)a 



J f V E W (tt 1 ^ £e^[i + e-«w]W) 

7Vs/ fe!,fe 2 \L rf {r}J m / 

x ^c^-fl + c-^-]^ 1 ^ ■ (— ^ ^e^"3[ 

\[ d {T}] n 2 J \[d{r}] n 3 

x [n FD (e fc2 ) - n FD (i/x„ 2 )] 



-i^„„l<i { r} 



7i FD (-£ fcl )n FD (-i^» 1 )n FD (i/i ?t3 ) + n FD (£ fcl )n FD (i/i„ 1 )n FD (-i^ K3 ) 
(iuj m e — £fej — i/i ni + i/x„ 3 )(£fc 1 — £fc 2 + i/x„ 2 — i/x n3 ) 



(E5) 



By examining the previous expression (|E5() one can sus- 
pect a technical difficulty in performing the summations 
over the imaginary chemical potentials (labelled as ni, 
n,2 and 713), due to their presence in the denominator. 
However a simplification occurs, because our study is re- 
stricted to the zero-temperature limit. Indeed the auxil- 
iary chemical potentials are proportional to the tempera- 
ture. Therefore their contribution to the denominator of 
Eq. I|E5(1 cancels out at zero-temperature. Of course such 
a simplification would not be valid if it was applied be- 



fore performing all the summations over the Matsubara 
frequencies: it would simply mean that the fermionic con- 
straint © is neglected, i.e., the considered Hilbcrt space 
includes unphysical fermionic states. However once the 
summations over the Matsubara frequencies have been 
computed, all the poles of the free propagators have been 
taken into account. Consequently the elimination of the 
imaginary chemical potentials from the denominator at 
this stage does not suppress any "hidden" property and 
is perfectly justified. We have to notice also in Eq. (|E5|) 
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that the contributions of the imaginary chemical poten- 
tials remain at zero-temperature in the exponential fac- 
tors, contrary to what happens in the denominator. In- 
deed the terms e 1/3/i ™ are purely trigonometric factors, 
independent of the temperature. It is probably feasible 
to extend the study to the finite temperature limit, for 
example by performing the summations over the imagi- 



nary chemical potentials numerically. This generalization 
is left for further studies. 

To obtain the expression of the physical self-energy, 
we perform the analytical continuation to the real axis: 
iu m e —> co + i0 + . Using the identities l|C9|l . (|C10|) and 
(|C1H in order to calculate the summations over the 
imaginary chemical potentials, we obtain at T = 



F, 



(3,0)a 



Ms 



CO — CO 



(co - e kl + iO+)(e fel - e k2 



riFD(efc 2 ) 



d {r} - 1 



d {r} - 1 



{r}J 



, (E6) 



r}J 



where co is substituted to iw TO e' after continuation to 
the real axis. We have now everything in hand to write 
the contribution to the electronic self-energy, at zero- 
temperature. We put together the spin factor (jDllfl and 
the frequency- momentum part i|E6() . and multiply them 
by two additional factors: (-1) due to the pseudofermion 
loop, and (iVj) the number of impurity spins. We get 



In the thermodynamic limit we adopt the standard cut- 
off schem o 26 ! 37 in energy for the conduction electron 
band: — D < e k < D, and perform the following sub- 
stitution 



-y 



— D 



de, 



^(3,0)a( ti; 7 w 



-1) • (Ni) • Spin 



J (3,0) £ 



x Ft 



(3,0)s 



where p is the density of states per spin and per channel, 
(E7) taken as a constant. We obtain finally 



(E8) 



(2S) 2 (2 -3N + N 2 ) +q 2 (2 + 3N + N 2 ) + 2S ■ g(4 - iV 2 ) +2S(-4 + 3N - 4N 2 + 2N 3 ) 



+q( - 4 - 3iV - AN 2 - 2N 3 ) + (2 + AN 2 + N 4 
x (Jk) 3 • TH ■ 5 a(T > ■ S(co - co ) 



(N + 2S-1)\ 2S 



{2S)l(q - l)\(N - q)\ AN 3 



x p~ I dsi / de 2 , , n , w r 

l-D J-d (w-ei+i0+)(ei-e 2 ) 



d in - 1 



d {r} - 1 



[d { 



r}J 



where m = Ni/Af s is the impurity density. 
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